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I— i' Context. Water vapour maser emission from late-type stars characterises them as asymptotic-giant-branch stars with oxygen-rich 

(~| , chemistry that are losing mass at a substantial rate. Further conclusions on the properties of the stars, however, are hampered by the 

^ |. strong variability of the emission. 

I ' Aims. We wish to understand the reasons for the strong variability of H2O masers in circumstellar shells of late-type stars. In this 

P ' paper we study RX Bootis and SV Pegasi as representatives of semiregular variable stars (SRVs). 

,y , Methods. We monitored RX Boo and SV Peg in the 22-GHz maser line of water vapour with single-dish telescopes. The monitoring 

^ , period covered two decades for RX Boo (1987 - 2007) and 12 years for SV Peg (1990 - 1995, 2000 - 2007). In addition, maps were 

, ^ , ' obtained of RX Boo with the Very Large Array over several years. 

Results. We find that most of the emission in the circumstellar shell of RX Boo is located in an incomplete ring with an inner radius 

I , of 91 mas (15 AU). A velocity gradient is found in a NW-SE direction. The maser region can be modelled as a shell with a thickness 

^ . of 22 AU, which is only partially filled. The gas crossing time is 16.5 years. The ring-like structure and the velocity gradient remained 

' stable for at least 1 1 years, while the maser line profiles varied strongly. This suggests that the spatial asymmetry is not accidental, so 

' that either the mass loss process or the maser excitation conditions in RX Boo are not spherically symmetric. The strong variability 

, of the maser spectral features is mainly due to incoherent intensity fluctuations of maser emission spots, which have lifetimes of the 

. order of 1 year. We found no correlation between the optical and the maser variability in either star. The variability properties of 

• ' the SV Peg masers do not differ substantially from those of RX Boo. There were fewer spectral features present, and the range of 

04 ' variations was narrower. The maser was active on the >10-Jy level only 1990 - 1992 and 2006/2007. At other times the maser was 

I either absent (<1 Jy) or barely detectable. 

00 , Conclusions. The variability of H2O masers in the SRVs RX Boo and SV Peg is due to the emergence and disappearance of maser 

■ clouds with lifetimes of ~1 year. The emission regions do not evenly fill the shell of RX Boo leading to asymmetry in the spatial 

il^ ' distribution, which persists at least an order of magnitude longer. 
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1. Introduction where mass-loss rates are high and they develop circumstellar 

dust and gas shells as AGB stars do. 

Water vapour maser emission of the 616 - 523 transition at , j- ^ 1 j- e u r\ • »u 

, „^ ^ , , N • ^ „ Immediately alter the discovery or H2O masers in the en- 

1 35-cm wavelength (22 GHz) is often found in the circumstellar ^j^^ j.^^ CMa, variability was noticed (Knowles 

shellsofoxygen-richstai-sontheasymptoticgiantbranch(AG ^1 ^ B^j^j p^^^ monitoring observations 

and in several red supergiants (RSGs). The AGB stars have in- ^^^^^ ^^^^ (Schwartz et al. [HB) evidence was 

termediate masses (M 1 - 8 Mq) and are usually long-period r j »u » »u u j » ■ a ■ 1 j • » 

. , , ^, ^ . , . /oT,, , s , found that the observed maser spectra varied in peak and inte- 

variables. They encompass semi-regular variables (SRVs), Mira » j a • i- j ui ■ i- * • u» 1 

• ,, .>x.T,Tr. r,,, , • .coT^x r gratcd flux, in linc profilcs and possibW in line-of-sight vcloc- 

variables and OH/IR stars. The classification of SRVs and Mira ^ ni /t-,- i- » 1 irm , t-u r » » j- • j 

. . , . , . , . , , ... . , , it\n (Dickinson et al. 1973). The earliest studies of long-period 

variables IS made in the optical regime and they are distinguished • ui * j »u ^ • u-i ,. n 1 » j 

,. , ^ , . . . ° . . , . , • , variables found that the maser variability was well correlated 

by the amplitude of their intensity variations. Due to their higher •., ... • .■ , j • • f j /c u » » 1 

^ , ^ ^,,„T^ • „ • , , . • , With variations in the optical and in the infrared (Schwartz et al. 

mass-loss rates OH/IR stars are optically faint and their variabil- mT/i, u »i » u \- u »»u 1 

, . r ./ • , , • 1974), but later observations showed that the regularity of the ra- 
tty IS observed in the infrared (see review by Habing 1996). .• • u ju . J /c 1 ; 1 ImoQl u 

\ J b dio variations had been overestimated (Engels et al. |1988l Cohen 

SRV, Mira and OH/IR stars form a sequence of increasing 
mass-loss rates, but it is not clear if this is also an evolutionary 

sequence or a sequence of increasing progenitor masses. RSGs ' I" this paper we use the expression 'line-of-sight velocity' instead 

are stars with masses M > 8 in a short evolutionary phase, °^ ''^'^'^^ ''^^°fy' ^^"^""'^ *e latter could be confused with the ve- 

locity radially from the star, i.e. the expansion velocity of the stellar 

envelope. All velocities in this paper are relative to the local standard of 

Send offprint requests to: A. Winnberg rest (LSR). 
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Table 1. The VLA observing schedule 



Date Conf. Sources 



lyyU reu. zo 


A 


KA Boo, U Her, (Jn 39. 






OH 83.4-0.9, VX Sgr 


1990 June 3 


A 


RX Boo, U Her, OH 39.7+1.5, 






OH 83.4-0.9, VX Sgr 


1991 Aug. 23 


A 


OH 39.7+1.5 


1991 Oct. 20 


BnA 


RX Boo, U Her 


1992 Jan. 10 


B 


VX Sgr, NML Cyg 


1992 Dec. 28 


A 


RX Boo, U Her, OH 39.7+1.5 


1995 Sept. 5 


A 


RXBoo 


1996 Dec. 10 


A 


OH 26.5+0.6 


1998 May 31 


A 


OH 26.5+0.6 



I1989I I. The first interferometric maps of the Mira variables R Aql 
and RR Aql taken at two epochs separated by 15 months also 
showed dramatic changes (Johnston et al. 1985). Knowledge of 
the properties of the variability and understanding of the reasons 
for it, however, are prerequisites for the interpretation of the ob- 
servations and for the derivation of basic stellar properties. 

Water vapour masers are potentially very sensitive tools to 
probe the physics and dynamics in parts of the circumstellar 
shells, whose radii range between ~5 and 50 AU in case of Mira 
variables and SRVs (Bowers et al. 119931 Bowers & Johnston 
[1994 ; Colomer et al.'2000', Bains et al.'2003', Imai et al. |2003] l. 
High-resolution mapping with interferometers shows that the in- 
tegrated maser emission is composed of a multitude of individ- 
ual clouds of sizes ~2 - 4 AU and a low filling factor of -0.01 
(Bains et al. l2003l l. In conjunction with the SiO masers located 
closer to the central star and with the OH masers in the outer 
parts of the envelope, H2O masers allow us to map the shell that 
is created by the strong mass loss (>10"^ M0yr"') on the up- 
per AGB. Information on the velocities and on the distribution 
of the maser spots can be used to reconstruct the geometry of the 
mass-loss process and sometimes even its history. 

Monitoring of H2O maser emission with single-dish radio 
telescopes, even covering more than a decade, provides only 
limited insight. The spectra mostly consist of many compo- 
nents, which often overlap in velocity. The profiles and inten- 
sities of H2O maser spectra therefore depend not only on tem- 
perature variations due to the pulsating star but also on other 
factors, including random intensity fluctuations of individual 
clouds (Engels et al. 119881 1. Early programmes followed only 
the strongest lines. Three years of monitoring led Cox & Parker 
( I1979I I to the conclusion that the masers are not stable for more 
than a few pulsation cycles. Later Little-Marenin and co-workers 
monitored a few Mira variables over a number of years and 
found complex variations in intensity with many components 
varying independently from one another (e.g. Little-Marenin et 
al. fT9#lD . Usuafly a phase lag of -0.3 (30 - 100 days) with re- 
spect to the visual maximum was observed (Little-Marenin et 
al. 119931 1. A Russian group monitored a number of stars with 
strong H2O masers for up to twenty years. They found that by 
and large the maser flux correlates with the visual light curve, 
and confirmed the presence of a phase lag. They attributed the 
maser variations and the phase lag to the passage of periodic 
shocks (Lekht et al. 1200 II and references therein). 

In order to improve the understanding of the properties of 
maser variability for different types of late-type stars, we started 
in 1987 a monitoring programme of several such stars using the 
Effelsberg and Medicina radio telescopes. The sample included 



SRVs, Mira variables, OH/IR stars and RSGs. In order to break 
the spatial degeneracy, which usually makes the interpretation 
of single-dish maser spectra difficult, we added interferometric 
observations using the Very Large Array (VLA). Of each object 
type we selected prototypical stars, which were observed during 
several epochs (see Table[T]). In this paper we present the results 
for two SRVs: RX Bootis (RX Boo) and S V Pegasi (S V Peg). In 
subsequent papers we will present results for the other types of 
late-type stars. 

RX Boo is a 9 - 11 mag SRV with a mean spectral type of 
M7.5 and a period varying between 340 and 400 days (Kukarkin 
et al. II97II 1. Olofsson et al. (2002) in their model fit involv- 
ing several CO transitions estimated a mass-loss rate of 6 x 
10"^ Mq yr"^ and an expansion velocity of 9.3 km s ' , whereas 
Teyssier et al. (120061 1. as a result of another model fit, obtained 
2 X 10"^ Mq yr~' and 7.5 kms"', respectively. The distance 
is 160+25 pc (The Hipparcos Catalogue; Perryman et al. |19971 l. 
Dyck et al. (119951 ). using an IR interferometer, measured the stel- 
lar diameter at 2.2 fim and found a uniform-disk (UD) diameter 
of 19±1 mas. They calculated an effective temperature of 3000 
+ 100 K. Chagnon et al. (120021) . observing in a band centred at 
3.8 fim, found a UD diameter of 21.0 + 0.3 mas. 

The H2O maser of RX Boo is one of the first masers of 
this kind detected in late-type stars (Schwartz & Barrett [T97 11 
Dickinson et al. 1973), and hence it is one of the best studied. 
Single-dish monitoring was performed by Schwartz et al. ( 1974]), 
Cox & Parker ( fT979] ). Berulis et al. ( iM3l ). and Hedden et al. 
( |1991i ) who all found a weak correlation of maser intensity with 
optical flux. Maps taken with the VLA show an overall extent of 
the maser shell of 150 - 230 mas (Bowers et al. |19931 Colomer 
et al. 2000), giving an outer shell radius of 12 - 18 AU or 9 - 14 
stellar radii. 

SV Peg also varies with ~2 mag amplitude (9-11 mag) 
and has a spectral type of M7 and a mean period of 144 days 
(Kukarkin et al.[1971j. The mass-loss rate is 3 x 10"^ Mq yr"' at 
an expansion velocity of 7.5 km s ' (Olofsson et al. 120021 ) and its 
distance is 205+45 pc (The Hipparcos Catalogue; Perryman et 
al. |1997i ). The H2O maser was discovered by Dickinson ( 119761 ) 
and was only occasionally observed before our monitoring ef- 
forts. No mapping has been done as yet. 

2. Observations 

Single-dish observations of the H2O maser line at 
22235.08 MHz were made with the Effelsberg 100-m and 
Medicin43 32-m telescopes at typical intervals of a few months. 
RX Boo was monitored over twenty years between 1987 and 
2007, and SV Peg in 1990 - 1995 and 2000 - 2007. VLAQ 
observations were made of RX Boo on four occasions in the 
period 1990 - 1992 and later in 1995 after a strong flare event 
(see Table[T]). This event will be discussed in a separate paper. 

2. 1 . Effelsberg observa tions 

With the 100-m telescope observations were made of RX Boo 
between 1990 February and 2002 June and of SV Peg between 

- The 100-m Effelsberg telescope is operated by the Max-Planck- 
Institut fiir Radioastronomie (MPIfR) at Bonn, Germany. 

^ The 32-m Medicina telescope is operated by the Istituto di 
Radioastronomia (IRA) at Bologna, Italy. 

"* The Very Large Array is operated by National Radio Astronomy 
Observatory (NRAO) from its subsidiary 'Array Operations Center' 
(AOC) at Socorro, New Mexico, USA. 
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1990 February and 1995 June. We used the 18 - 26-GHz re- 
ceiver to observe the 6i6 — » 523 transition of the water molecule. 
The receiver had a cooled maser as pre-amplifier and passed 
left-hand circularly polarised radiation. As circumstellar wa- 
ter masers are found to be unpolarised to limits of a few per- 
cent (Barvainis & Deguchi 1989), significant flux losses due to 
the polarization properties of the receiver are not expected. At 
1.3-cm wavelength the beamwidth is ~40" (FWHM). We ob- 
served in total power mode integrating ON and OFF the source 
for 5 min each. 'ON-source' the telescope was positioned on 
RX Boo at the coordinates a = 14'>24™ 1P.6, S = +25° AT 13" 
(J2000) and on SV Peg at a = 22^ 05'"42U, 5 = +35° 20' 54" 
(J2000), while the 'OFF-source' position was displaced 3' to the 
east of the source. Typical system temperatures were 60 - 120 K 
depending on weather conditions and elevation. The typical rms 
noise-level in the spectra is -0.2 Jy. 

The backend consisted of a 1024-channel autocorrelator 
Observations were made with a bandwidth of 6.25 MHz, cen- 
tred on the systemic line-of-sight velocity of the star. The ve- 
locity coverage was ~80 kms"' and the velocity resolution 
0.08 km s-'. 

Standard procedures were used to reduce the spectra. First a 
baseline was fitted using the channels not affected by emission, 
then the calibration was applied. Finally, spectra from several 
integrations of the same source were averaged, if they were taken 
during the same day. Hanning smoothing to 0.16 kms ' was 
applied only to noisy spectra. 

The flux calibrations were made against NGC 7027 (5.9 Jy) 
and Mon R2 (6 Jy), and corrections for atmospheric attenuation 
and elevation were applied using a standard gain curve deter- 
mined under good weather conditions. Deviations from the stan- 
dard gain curve were less than 10%, except at low elevations 
where deviations up to 30% were measured. Losses due to point- 
ing are at most secondary effects, as pointing checks and correc- 
tions were regularly carried out on strong continuum sources, 
showing deviations from the nominal position of at most 10". 

2.2. Medicina observations 

Observations at Medicina were carried out between 1987 March 
and 2007 January on RX Boo and between 1990 February and 
1995 January and again between 2000 October and 2007 January 
on SV Peg. The Medicina 32-m telescope (HPBW at 22 GHz 
~ 1.'9) is set up primarily to be used for VLBI measurements, 
and therefore the front-end is typical for such work. For line ob- 
servations at 22 GHz only the right-hand circular polarization 
output from the receiver is available. 

Over the years the system has undergone various improve- 
ments. In 1989 realignement of the antenna surface resulted in 
an improvement of the efliciency at 22 GHz from 17 to 38%. 
The backend was a 512-channel digital autocorrelator; in 1991 
the number of channels was increased to 1024. In the same year a 
HEMT amplifier replaced the GaAs FET front-end, reducing the 
system temperature (120 K at zenith under good weather condi- 
tions). Also in 1991 an active sub-reflector control increased the 
gain at low and high elevations. The available bandwidth varied 
between 3.125 MHz and 25 MHz. Up to 1997 the spectral res- 
olutions used corresponded to 0.082 kms"' or 0.165 kms"'. In 
1997 the VLA-1 chip correlator was replaced by one based on 
the NFRA correlator chips and boards (Bos |199Tl ); in the stan- 
dard configuration the latter has 2048 channels and a maximum 
bandwidth of 160 MHz. However, to maintain consistency with 
the older observations we have only used a 10 MHz band and 
1024 channels, giving a velocity resolution of 0. 132 km s"' . 
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Table 2. VLA map specifications for RX Boo 



UaiC 




LLrD vv 




rms 


<^ IM 
J / 1\ 




maj.a. 


min.a. 


p. a. 








(") 


(") 


D 


(Jy/b.) 




1990 Feb. 


0.090 


0.077 


86.23 


0.015 


3700 


1990 June 


0.082 


0.078 


40.12 


0.014 


1700 


1991 Oct. 


0.305 


0.092 


-76.62 


0.025 


190 


1992 Dec. 


0.077 


0.072 


-23.67 


0.027 


130 


1995 Sep. 


0.127 


0.084 


-67.27 


0.009 


580 



maj.a.: half-power beam width (HPBW) for the major axis of the best- 
fit three-dimensional Gaussian component to the synthesised beam 
min.a.: HPBW for the minor axis 
p.a.: position angle of the major axis (E of N) 

rms: the root-mean-square noise fluctuations in signal-free channels in 
units of Jansky per beam area 

S/N: signal-to-noise ratio or 'dynamic range' in the channel with 
strongest signal 



The telescope pointing model is typically updated a few 
times per year, and is quickly checked every few weeks by ob- 
serving strong maser sources (e.g. W3 OH, Orion-KL, W49 N, 
and/or Sgr B2 and W51). The pointing accuracy was always bet- 
ter than 25"; from about mid 2004 the rms residuals from the 
pointing model were of the order of 8"- 10". The antenna gain 
as a function of elevation was determined by daily observations 
of the continuum source DR 21 at a range of elevations; we as- 
sume a flux density of 16.4 Jy (after scaling the value of 17.04 Jy 
given by Ott et al. ( 1994) for the ratio of the source size to the 
Medicina bearr0). The daily gain curve was determined by fitting 
a polynomial to the data, which was then used to convert antenna 
temperature to flux density for all spectra taken that day. From 
the dispersion of the single measurements around the curve we 
find the typical calibration uncertainty to be 19%. On the few 
days that no separate gain curve was measured, we have applied 
the one closest in date, and estimate a corresponding calibration 
uncertainty of 7%. The average overall calibration uncertainty is 
therefore estimated at ~20%. 

Observations were taken in total power mode, with both ON 
and OFF scans of 5-min duration. The OFF position was taken 
1° W of the source position. Typically, 2 ON/OFF pairs were 
taken. The typical rms noise-level in the spectra is ~1.5 Jy. 
Finally, no observations were taken between 1996 April and 
1997 February, due to extensive maintenance work on the tele- 
scope. 

For the purpose of studying the long-term variability of the 
maser emission, spectra of the same source taken less than 4 days 
apart were averaged. During the period that the sources were also 
monitored at Effelsberg it may happen that we have spectra from 
both telescopes within 4 days. In these cases we always used the 
EflFelsberg spectra because of their much lower noise level. With 
very few exceptions both spectra showed the same line profile 
and flux densities (within 10%) of the individual emission com- 
ponents, providing an a posteriori check on the correctness of 
the calibration procedures of both telescopes. 
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2.3. VLA observations 

RX Boo was observed with the VLA on four occasions between 

1990 February and 1992 December (Table[T]|. An additional ob- 
servation was made in 1995 September after the occurrence of a 
strong flare, which will be discussed in a separate paper. All 27 
antennas were used yielding synthesised beamwidths as speci- 
fied in Table|2] These beamwidths depend mainly on the config- 
uration, i.e. the extent of the array. For three of the four epochs 
we used the largest extent (A; see Table [T]), while the October 

1991 observations were carried out with a hybrid configuration 
(BnA). 

We chose a receiver bandwidth of 3.125 MHz to obtain a 
total velocity range of 42 kms"' and the bandwidth was split 
into 64 channels, yielding a velocity resolution of 0.66 kms"'. 
Data from the right and left circular polarization modes were 
averaged. Typical integration times were 30 min on the star and 
12 min on the phase calibrator J1407 -1-284 with a sampling time 
of 30 sec. Flux calibration was obtained relative to 3C286 that 
was assumed to have a flux density of 2.55 Jy and 3C84 was 
used to correct for the bandpass shapes. 

The data were reduced with AIPS (Astronomical Image 
Processing SysterrQ) using the spectral line data routines. To 
improve the sensitivity as well as the dynamic range, self- 
calibration was applied. Care was taken to select an unresolved 
feature as phase reference. It turned out that bright features of- 
ten are extended and we made several trials to find the best un- 
resolved feature. Usually we used a point source model for the 
first iteration. The final sensitivity reached on line-free channels 
was 9-27 mJy/beam depending on weather conditions and in- 
tegration time (see Table|2]i. 

3. Single dish data of RX Boo 

3.1. Definitions 

To interpret the single dish maser profiles we will use several 
terms, whose definitions are given here. The integrated maser 
emission from a circumstellar shell is a superposition of emis- 
sions arising from many physical entities {maser clouds). These 
emissions are in the form of narrow maser lines, which emit at 
the velocity of the cloud projected on the line-of-sight. Observed 
with a single radio telescope all the maser lines together form a 
maser spectrum with usually several peaks. The shape of the of- 
ten contiguous spectrum is named maser profile. The data anal- 
ysis aims at a decomposition of the spectrum into individual 
maser lines. However, due to blending in velocity space, the 
analysis rarely results in individual lines but in maser features, 
which are composed of two or more maser Unes. These maser 
features might be traced over some period of time in several 
spectra and may reappear at later times perhaps with a slightly 
different velocity (^ 0.5 kms"'). They might be identified with 
real maser clouds persisting over this period of time. However, 
most likely the features are composed of several real clouds 
with velocities close to one another but spatially separated. We 
will use the interferometric maps to identify these (groups of) 
maser clouds, which we name more generally maser spatial 
components. For accounting purposes all spectral features were 
grouped according to their velocities into maser spectral com- 
ponents. At different times within the monitoring period of 20 

^ In previous works and up to 2003 April a flux density of 18.8 Jy 
was adopted for DR 21, following Dent (19721. In this paper, all data 
have been recalibrated using the value given by Ott et al. l fT^94l l. 

* http://www.aoc.nrao.edu/aips/ 



years different spectral features contribute to a spectral compo- 
nent, reflecting the asynchronous formation and dissolution of 
the contributing maser clouds. 

3.2. H2O maser spectra 

Our monitoring of RX Boo covers almost 21 years, from 1987 
March to 2007 January with several (typically 4-5) obser- 
vations per year. Depending on telescope (i.e. Effelsberg or 
Medicina), date of observation and integration time the rms sen- 
sitivity of the observations was very inhomogeneous ranging 
from 0.1 to 6 Jy. All spectra taken are shown in the Appendix 
(Fig. Al). Sample spectra showing typical profiles are given in 
Fig.[T] A general view on the properties of the masers and their 
variations is given in Fig.|2] Here the spectra are arranged along 
a time axis and the flux densities in between observations are in- 
terpolated to obtain an approximated continuous variability pro- 
file. 

The H2O masers of RX Boo covered a constant velocity in- 
terval between -5 and +5 kms"', although at particular veloci- 
ties the maser emission could drop below the detection threshold 
for extended periods of time (months to even years). In the time 
range 1987 - 1989 (t ^ 1000 in Fig. ^ RX Boo displayed a 
strong maser feature at 2.5 kms"' reaching flux density levels 
of a few hundred Janskys. In 1995 (f ~ 2900) a flare was ob- 
served at a velocity of 3.4 kms"' with a peak flux density of 
>1500 Jy. During the rest of the time typical flux density lev- 
els were <100 Jy. The maser never completely disappeared, al- 
though occasionally the flux density of the strongest emission 
peak was only a few Janskys. 

Although strong flux variations were observed, leading to a 
great variety of spectral profiles, they did not occur at random 
velocities within the observed velocity interval. With a few ex- 
ceptions the emission at vios > km s"' (observed line-of-sight 
velocity) was always stronger than that at vios < kms"', with 
the strongest peak occurring between 1 .5 and 3.5 km s ' . 

3.3. Line profile analysis 

The H2O maser spectra of RX Boo were decomposed into sepa- 
rate features by fitting multiple Gaussian line profiles. We gener- 
ally used an iterative procedure in which we first fitted and sub- 
tracted the strongest features in each spectrum and then extracted 
the weaker features from the residual spectrum. The number of 
extracted features depended on the quality of the spectra as well 
as the number of distinct peaks. Blending of maser lines is a 
major obstacle and the extracted features may consist of two or 
more lines. Depending on the strength of the contributing lines 
the central velocity of the extracted feature may vary. As a rule 
of thumb, the FWHM of strong features (visible as distinct peaks 
in the spectra) is ~1 kms"', and therefore features with FWHM 
2 kms"' are most probably blends. Comparison with adja- 
cent spectra often shows that blended features split when the 
relative intensities of the contributing lines change. In complex 
situations when Gaussian fitting failed, features were identified 
by eye and velocities and flux densities measured by hand. We 
assumed that maser features in adjacent spectra with velocity 
differences 50.5 kms"' belong to a unique spatial maser com- 
ponent in the shell of RX Boo, which varied in intensity with 
time. All features were accordingly assigned to 1 1 maser spec- 
tral components labeled A, B,..., K. These components are listed 
in Tables lA.ll and IA.2I in the Appendix. For each spectrum we 
give the Gregorian and Julian dates, its rms noise level in Jansky, 
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-10 10 -10 10 -10 



31-JAN-1989 673 09-MAR-1995 2901 ia-NOV-2003 6077 




-10 10 -10 10 -10 10 

22-DEC-1992 2094 06-MAY-1997 3690 14-FEB-2006 6896 




-10 10 -10 10 -10 10 



V|os (km/s) 

Fig. 1. Sample H2O maser spectra of RX Boo showing typical profiles starting in 1987 and ending in early 2007. Note the varying 
flux density scale. 
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RX BOO 




-5 5 



Vtos (km/s) 

Fig. 2. Flux density of the H2O masers of RX Boo as a func- 
tion of velocity and time. The time axis covers the years 1987 - 
2007. The time stamp t = corresponds to 1987 March 30 (J.D. 
2446885). The dotted horizontal lines mark the actual observa- 
tions averaged in 4-day intervals. Flux densities in between the 
observations were interpolated linearly. The spectra were resam- 
pled to a resolution of 0.3 km s ' . The dashed vertical line marks 
the stellar systemic velocity (1 kms"') as given by Teyssier et 
al. (2006). Contours are drawn at levels 3.00, 8.49, 24.03, 68.02, 
192.51, 544.85 Jy (logarithmic scale). These six levels are indi- 
cated by black and white lines in the flux scale along the right 
side of the map. (For a colour version of this figure, see the elec- 
tronic edition.) 

the Si in 10"^^ W m"^, and the velocities (Vp) and peak flux den- 
sities (Sp) for each component. Occasionally, features were con- 
sidered as blends and were assigned to two or more components. 
Values marked by ":" and upper limits were measured with the 
cursor on the computer screen. 

3.4. Line fitting results 

Tables lA.ll and IA.2I show that the maser features at 
vios < kms"' were weaker than those at vios > kms ' 
over the full monitoring period. One spectral component, C, at 
-2.7 kms"', was detected at almost every epoch, whilst A, B, 
D and E were only detected sporadically. This may be largely 
a sensitivity effect, because many spectra did not have a suffi- 
cient signal-to-noise ratio to allow detection of the fainter fea- 
tures on the level of a few Janskys. In the red part of the maser 
profile (vios > kms ') heavier blending of lines made the 
unambiguous assignment of features to the components F - K 
a non-trivial task (Table IA.2b . The features assigned to com- 
ponent I, for example, show a systematic velocity shift with 
Vp= H-2.7+0.3kms-' until 1992 and -1-3.2+0.3 km s"' after 1994. 
Most likely the emission contributing to component I came from 
(at least) two maser clouds, of which one of them dominated 



the emission during the first time range and the other one during 
the second time range. Similar shifts are observed among the 
other spectral components. The direction of the shifts does not 
show any systematics so that we regard these shifts generally 
as a result of blending by two or more maser lines contributing 
to a spectral component. Some components show a rather large 
scatter of velocities, as for example component G with 0.4 < 
Vp< 1.9 kms"'. These "shifts" are almost random and are also 
not well explained as movements of a physical maser cloud, but 
rather as a result of blending of several maser lines in the ve- 
locity interval 0-2 kms"' and varying relative intensities. The 
spectral components F and G show therefore some overlap in 
velocity, indicating that they are not well defined. Spatial reso- 
lution of the maser emission is required to overcome the limited 
interpretational power of the single dish maser profiles. 

However, it is obvious from Fig. |2]as well as from the run 
of velocities in individual spectral components (Tables [ATI and 
IA.2I 1 that the lifetimes of individual components are limited on 
the scale of months up to a few years. The lifetimes will be ad- 
dressed in more detail in Sect. 14.51 

4. Interferometric data of RX Boo 

Interferometric observations of the maser emission aim at map- 
ping the maser clouds and breaking the spatial degeneracy. Due 
to the limited spatial resolution of the instrument this may be 
only partially successful, and the data analysis of the images will 
yield maser spatial components, which will be superpositions of 
several maser clouds close to each other in space as well as in 
velocity. 

4. 1 . Cliannel-map analysis 

The VLA interferometric data consist of one data cube for each 
of the four observations, containing 63 maps each. The maps are 
separated in velocity by 0.658 kms"'. An example is given in 
Fig. [3] where maps of the 22 channels with detected signal are 
displayed from the first VLA epoch (1990 February). Except for 
the central velocity range the maser emission consists of point- 
like and unresolved spots. In order to single out maser compo- 
nents the data cubes were analyzed within AIPS in a three step 
process. First, the individual channel maps were analyzed one by 
one, then components were identified by comparing neighbour- 
ing channel maps, and finally these components were verified in 
velocity space. 

The analysis of the channel maps was made by fitting up to 
four two-dimensional Gaussian components with elliptical cross 
sections to the emission spots. The results of this analysis tend 
to be insecure as soon as spots are blended spatially as well as 
in velocity. The positions of the fitted Gaussian components are 
obtained with an accuracy of ~0.02 pixel (0.4 mas) for strong 
emission spots (>1 Jy). Comparing the positions of such spots 
extending over several channels yields deviations <0.2 pixel 
(<4 mas) between the positions of the same spot in different 
channel maps. The positions of fainter emission spots have er- 
rors up to 0.5 pixel (0."01). Very faint spots (<300 mJy) can be 
localised by summing over several channels, usually three. This 
is guided by the experience that spatial components extend in 
velocity by not more than ~3 channels (~2 kms"'). Often the 
Gaussian components fitted to weak spots from merged channel 
maps tend to be extended, indicating that these spots are actu- 
ally blends of several weak emission spots. Due to high signal- 
to-noise ratios the positional accuracy obtained from Gaussian 
fits is much better than the resolution of the interferometer of 
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Fig. 3. Sample H2O maser images of RX Boo from 1990 February. Each image is separated from its neighbours by the line-of-sight 
velocity interval 0.658 kms"^. The synthesised FWHM beam size (major axis: 0."09; minor axis: 0."08; position angle of the major 
axis: 86°) is shown in tlie uppermost left image (7.9 kms"'). The images are oriented along right ascension and declination and 
the angular scales are relative to the position a = 14*^24" 11^6, 6 = -i-25°42' 14" (J2000). Brightness contours are -0.05 (dashed 
contour), 0.05, 0.5, 5, 50 Jy per beam area (1.45 10"^^ sr). 
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~4 pixel (FWHM ~ 0."08), allowing us to determine the relative 
spatial separation of spots down to 4 - 10 mas. 

4.2. Spatial component identification 

Subsequently, neighbouring channel maps were compared and 
spatial maser components were identified by looking for spa- 
tially coincident Gaussian components. Spatial maser compo- 
nents were accepted as real, if they had Gaussian components in 
at least two neighbouring channel maps. Gaussian components 
present in only one map are always weak, and are not consid- 
ered further. They might be real maser components but cannot 
be verified. At last we made sure that the list of identified spa- 
tial maser components can be identified also in velocity space. 
For this purpose 'slices' along the velocity axis through the cube 
were obtained at the spatial positions of the components. Usually 
the components stand out as distinct peaks in these 'slices', con- 
firming their reality. Spatial components blended in velocity by 
another strong component, reveal themselves by adding consid- 
erable asymmetry to the velocity profile of the major compo- 
nent. Asymmetric velocity profiles of components is a strong 
indication of blending and may reveal new components that are 
separated marginally in velocity but not spatially within the res- 
olution of the maps. The final list of components consists of all 
components verified spatially as well as in velocity 

Typically 10 to 15 spatial maser components were identified 
in each VLA observing epoch. These components are listed in 
Table [3] For the epoch 1990 June the components can be com- 
pared to those found by Colomer et al. 62000 ). who used the 
VLA to observe RX Boo one day apart from our observation. To 
determine the spatial components they used a 3-D Gaussian fit- 
ting programme developed to analyze directly the spectral data 
cubes. They find about the same number of components (15 vs. 
14), and their Gaussian-fitted intensities and the overall spatial 
distribution of their components are compatible with our results. 
The rms scatter in intensities and positions for components in 
common are 2 Jy and 30 mas, respectively. The positional de- 
viations are a factor >3 higher than the relative positional ac- 
curary of 4 - 10 mas derived above. A possible cause for the 
scatter is the extensions of the maser clouds. Bains et al. ( I2003I I 
found sizes of 2 - 4 AU for the Mira stars U Ori and U Her 
Adopting these sizes for RX Boo, angular sizes of ~40 mas 
are predicted, much larger than the implied positional accuracies 
from the Gaussian fits. We suspect that the two fitting methods 
treat the eff'ect of blending of the extended maser clouds differ- 
ently, spatially as well as kinematically. 

To search for continuum emission of the star we formed the 
average of all channels without line emission, yielding rms fluc- 
tuations of the order of 5 mJy/beam. The non-detection is in 
agreement with the expected blackbody flux density of ~0.3 mJy 
from a star with a diameter of 19 mas and a photospheric tem- 
perature of 3000 K. 

4.3. Alignment of ttie maps 

In order to identify spatial components present over two or more 
observing epochs we aligned the four maps to a common ori- 
gin. The maps of 1990 were very similar showing the majority 
of components to be present in both maps and making the align- 
ment trivial. Several spatial components in common were found 
also in the 1991 and 1992 maps, however, the components in 
common with the 1990 maps were scarce. Finally, we chose the 
strong spatial components Gl, HI and 11/12 (see Table |3]l to ob- 



tain an alignment of aU four maps. After alignment these com- 
ponents scattered in position by ~10 mas. 

Spatial coincidences among the weaker components were 
searched for in the aligned maps. After subtraction of coincident 
components we ended up with 25 different spatial components 
of which 64% were present in at least two maps. Positional de- 
viations between maps were typically of the order of 30 mas, 
although in a few ambiguous cases we accepted as coincidences 
also components with deviations up to 60 mas (cf. El, H3, Kl 
in TableO. The velocity and spatial resolutions were clearly in- 
sufficient to avoid such ambiguities. 

According to their velocity, the merged spatial components 
were assigned to the eleven spectral components A...K. The 
merged components are listed in Table[3] which gives the spatial 
component (Spat. C), the VLA observing epoch, the line-of- 
sight velocity and peak flux of the spatial component identified, 
the spatial offsets (Xoff and Yoff) from the adopted map cen- 
tre (see below) and the associated spectral component (Spec. C.) 
from Tables|AT|and|Aj2]in the Appendix. 

4.4. Distribution of spatial components 

The distribution of spatial components from all maps after align- 
ment is shown in Fig.|4l The flux density of the spatial compo- 
nents is represented by the sizes of circles centred on the position 
of the components. Due to the large range in flux (see the third 
column of TableO we have chosen a logarithmic relation: 

£/ = 20 log 5 4-26 (1) 

where d is the diameter of the circle in milliarcseconds and S is 
the flux density in Janskys. 

In this map the most prominent feature is the incomplete 
ring-like structure where most of the spatial components reside. 
The south-western part of the ring contained no emission dur- 
ing all four observing epochs covering 3 years. The second ob- 
vious feature is a velocity gradient with most of the redshifted 
emission (v'los > 2 kms"') located in the south-eastern part and 
most of the blueshifted emission (vios < -1 kms"') in the north- 
western part of the shell. Although containing fewer spatial com- 
ponents, also the VLA map taken in 1995 September shows the 
velocity gradient and is compatible with the incomplete ring ge- 
ometry. 

We interpret the ring-like structure as outlining the inner 
boundary of tangential masers. Therefore, the position of the 
central star should be close to the centre of the circle that fits 
best to the components in this ring-like structure. Consequently, 
we have made a least-squares fit of a circle to all components 
except A2, F2, HI and 13 (of all epochs), which are situated off 
the ring-like structure at about 1 1 o'clock and 1 o'clock. The ra- 
dius of this circle is 91 +6 mas and the statistical errors of the 
centre position are of the same order of magnitude (6-7 mas). 
We have then moved the origin of the distribution from an earlier 
arbitrary position to this centre. The spatial offsets Xoff and Yoff 
listed in Table [3] are given relative to this centre. The choice of 
the stellar position is discussed further in Sect. 18.31 

4.5. Cross correlation of single-dish and interferometric data 
1990-1992 

As suspected already in Sect. 13.41 the spectral components are 
actually blends of two or more spatial components. During 1990 
- 1992 many spectral components were made up of coeval spa- 
tial components (e.g. Al and A2 in Table [3]). At later times the 
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Fig. 4. All the maser components listed in Table [3]plotted on the sky. Each component is represented by a symbol and a circle whose 
diameter is logarithmically proportional to the flux density of the component (see text for the mathematical expression). The epochs 
are represented by different symbols: 1990 February by small circles; 1990 June by asterisks; 1991 October by plus signs; 1992 
December by crosses. The circles around the components are colour coded according to the line-of-sight velocity of the component 
(vios in the text). (For a colour version of this figure, see the electronic edition.) A circle (dashed) has been fitted to the components 
on the arc (see text for a detailed list of components) and the origin has been moved to the centre of this circle (TableO. The filled 
dark red circle at the origin symbolises the central star with a diameter of 19 mas. 



same spectral component may come from emission of entirely 
different spatial components. The drift in velocity seen in com- 
ponent C, for example, is accompanied by a change of the emis- 
sion region within the shell (see spatial components CI - C4 in 
Table[3]l. Some components, as for example the most redshifted 
emission (spectral components J and K), are only temporarily 
present. 

While the spatial components detected in the 1990 February 
map were mostly redetected in 1990 June, many of them had dis- 
appeared in 1991 October Also there were fewer coincidences 
of spatial components between the 1991 and 1992 maps, com- 
pared to those between the 1990 maps. This indicates that the 
typical lifetime of individual maser clouds is in the range 0.5 - 1 



year. However, spectral components J and K could be detected in 
the single-dish spectra over 1 .3 and 2.5 years, respectively (Table 
IA.2I 1. We conclude therefore that there is a range of lifetimes of 
individual maser clouds of 0.5 - 3 years. The regions within the 
shell able to sustain maser emission at least from time to time, 
seem to have a longer lifetime, otherwise the stability of the dis- 
tribution of maser components over ~3 years cannot be under- 
stood. As the maser components we detected with the VLA are 
probably made up of several maser clouds, the quoted lifetimes 
are certainly upper limits for the lifetimes of single maser clouds. 
Indeed VLBA maps taken in 1995 by Mai-vel (1996) show a 
much larger number of spatial components, which would have 
been impossible to resolve with the VLA. 
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4.6. Proper motion 

Over the period of the VLA observations (2.8 years) an indi- 
vidual cloud with the maximum possible velocity perpendicular 
to the line of sight, i.e. 8.4 kms"' (see Sect. 15.21 ). would have 
moved ~30 mas. With much smaller line-of-sight velocities and 
lifetimes mostly considerably shorter than 3 years no detectable 
proper motions of spots within the H2O maser shell were ex- 
pected. However, the proper motion of the whole system was 
likely to be detected. 

To test for the presence of proper motion of RX Boo be- 
tween the 1990 and 1992 VLA observations we searched for 
spatial maser components present on both occasions in maps 
obtained from the calibrated dataset, in which phase calibration 
was made with the use of the phase calibrator J1407 4-284, i.e. 
before self calibration was applied. Only spatial components Gl 
(0.2 - 0.8 km s-i) and HI (1.4 - 2.1 km s-i) could be identified 
unambiguously in the 1990 and 1992 maps, while they could 
not be recovered in the 1991 October maps, due to decreased 
resolution (BnA configuration) and lower signal-to-noise ratio. 
To measure the positions of the components we used the AlPS 
task IMFIT. Positional differences between the epochs were av- 
eraged for both components yielding Aa - +49 ±11 mas and 
Ai5 - -159+6 mas. These differences are attributed to a common 
movement of maser sites and central star. The proper motion of 
RX Boo was then Aa - -1-18+4 mas/yr and AS = -59+2 mas/yr 
For comparison, the proper motion determined by Hipparcos is 
Aa = +21.74 + 0.75 mas/yr and A6 = -49.70 + 0.69 mas/yr 
(Perry man et al. 1997). Using the geometry suggested in Fig.|4] 
we determined the likely position of RX Boo for 1990 February: 
a = 14''24" 1F.52 + 0501, 6 = +25° 42' 13'.'7 + O'.'l (J2000). 



5. The circumstellar envelope of RX Boo 

5. 1 . The projected structure 

In an isotropically expanding envelope the relationship between 
the line-of-sight velocity and the projected distance from the star 
is an elliptical one: 



(2) 



where r is the distance, is the projected distance, v is the ex- 
pansion velocity, vios is the line-of-sight velocity and v, is the 
line-of-sight velocity of the star (+1 kms"'). Therefore, when 
plotting (vios - Vt) against rp for spatial components in such an 
envelope where, in addition, all the masers are situated in a thin 
spherical shell (like, for example, OH masers in an OH/IR star), 
the points will be close to (half) an ellipse with its centre at the 
origin and its axes coinciding with the coordinate axes. However, 
if the components are spread out within the envelope, as is ob- 
viously the case for the H2O masers in RX Boo, the plot will 
become a scatter diagram without any discernable ellipse. On 
the other hand if, for example, there is a minimum distance from 
the star where masers can exist and the component density is 
high, there will be a void in the lower part of the diagram with a 
boundary that will fit well to an ellipse. 

In order to test this method of analysis on RX Boo we 
plotted the absolute value of the line-of-sight velocity of the 
spatial components in Table [3] versus the projected distance 

( VXoff^ + Yoff"^) from the star (Fig.|5j. By choosing |vios - v»|, 
the density of points in the plot is doubled and the chance of 
seeing a possible boundary is increased. The fact that no masers 
are seen close to the star in Fig. |4] and no spatial components 
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Fig. 5. The absolute value of the line-of-sight velocity relative 
to that of the star (1 kms"') plotted against the projected dis- 
tance of all the spatial components listed in Table [3] from the 
origin (the central star). The epochs are represented by the same 
symbols as in Fig. |4] Five of the components (marked by sur- 
rounding circles) fit accurately to an ellipse (dotted line) centred 
on the origin and with its axes along the coordinate axes of the 
diagram. If the expansion of the stellar envelope is isotropic, this 
suggests that there exists an inner distance (14-16 AU) at a cer- 
tain expansion velociy (4-6 km s~') where water vapour masers 
first appear in the expanding gas. 



are found at projected distances smaller than ~9 AU in Fig. |5] 
are solid proof that there are no masers radiating in the radial 
direction ('radial masers') in the circumstellar envelope (CSE) 
of RX Boo. Most of the masers in fact are radiating in direc- 
tions that probably are close to tangential ('tangential masers'). 
They can be seen forming a small cluster around about 15 AU 
(94 mas) from the star at velocities close to zero in Fig. |5] cor- 
responding to several of the components along the dashed circle 
in Fig. |4] with radius of 91 mas. However, there is a mixture of 
line-of-sight velocities among the components close to this cir- 
cle, which is clearly seen in Fig.|5] namely the roughly vertical 
extension of the above mentioned cluster towards higher (in an 
absolute sense) line-of-sight velocities. There are about 10 com- 
ponents spread out at larger distances (rp > 20 AU) from the star 
visible in both figures. They must be located in the gas that has 
been carried out to these distances in the outflow. The fact that 
the line-of-sight velocities are not very high suggests that these 
masers too are more or less tangential masers. 

In Fig. |5]a quarter of an ellipse has been drawn. This repre- 
sents a least-squares fit to five of the data points in the diagram, 
marked by circles around the symbols. The best-fit values of its 
axes are 14.6 +1.4 AU and 4.8 ± 0.9 kms"'. We interpret this 
well-fitting ellipse as an indication of the existence of a mini- 
mum distance for water vapour masers of 14 - 16 AU in the CSE 
of RX Boo. At this distance the expansion velocity probably is 
in the range 4-6 km s"' . 



5.2. The 3-D structure 

By assuming a plausible expansion law it is possible to get a 
rough idea of the three-dimensional structure of the envelope. 
Supported by theoretical models of radiation-pressure-driven 
mass loss from late-type stars (cf. Goldreich & Scoville 119761 
and Bowen [T988I I we assume that the expansion law is expo- 
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Fig. 6. The model expansion velocity as a function of the dis- 
tance from the star (Eq. O. It approaches the terminal velocity 
(8.4 kms"') asymptotically (horizontal straight line). The dots 
along the expansion law curve symbolise the distances and ve- 
locities of the maser components. 

nential in nature and leading asymptotically to the final velocity, 



V = Vf {1 - exp [-k{r - ro)]} 



(3) 



where A: is a scaling factor and ro is a radial shift. 

Eq. (O thus contains three constant parameters (vf , k and ro) 
and in order to find plausible values for them we need to estimate 
the coordinate values for three points on the expansion law. The 
value of Vf (the expansion velocity at r = oo) can be estimated 
from (sub)-millimetre CO data, because it is believed that CO 
emission originates from gas throughout the CSE to very large 
distances from the star. We have two estimates of the terminal 
expansion velocity available: 9.3 kms ' (Olofsson et al. |2002| l 
and 7.5 kms"' (Teyssier et al. 2006). Both these values are the 
results of complicated model fits involving several rotational CO 
lines. Therefore we are unable to judge the correctness of these 
results and we simply assume that they reflect the uncertainty 
inherent in the two models that give rise to them. Therefore we 
adopt the mean value of 8.4 kms ' for the value of the final 
expansion velocity. 

In addition, we have a rough idea of the expansion velocity 
at one distance from the star (4.8 kms"' at 14.6 AU) from the 
analysis of the diagram in Fig. |5] 

Finally, we assume that the expansion law is purely expo- 
nential all the way down to the photosphere, which is stationary, 
and where the expansion velocity is zero. 

Thus the values of the parameters k and ro can be deter- 
mined from the condition that the stellar wind passes the po- 
sitions (ri,vi) = (14.6,4.8) and (r2,V2) = (1.5,0): 



In Vf - In (vf - vi) 
r\ - ro 



and 



ro 



(4) 



(5) 



resulting in A: = 0.065 AU" and ro = 1.5 AU. 

Now, the observables are the line-of-sight velocity relative 
to the star (vios - v,) and the projected radial distance rp. Their 
relations to the physical expansion velocity v and radial distance 
r are: 



rr, = r cos ( 



(7) 



V, = V sin 6 



(6) 



where 6 is the angle between the normal to the line of sight and 
the radius from the star to the maser. Therefore, for every maser 
component the angle can be determined from the equations (|3]l, 
© and 0. 

In Fig. |6] the expansion law is shown graphically together 
with the maser component values determined by the model and 
represented by filled circles. They are shown mainly as an indi- 
cation of the distance range for the maser components (~14 to 
-36 AU). 

The 3D-distribution of the maser components and their asso- 
ciated velocities is shown in Fig. Q as seen from three different 
directions: on the sky (a; similar to Fig.|4]but with the projected 
expansion velocities added), 'from the side' (b) and 'from above' 
(c). The dominance of tangential masers is seen clearly in Fig.|2] 
(b) and (c). The obliquity of the 'maser arc' seen in Fig.|2](c) is in 
accordance with the preference of negative line-of-sight veloci- 
ties on the lower-RA side of the star and of positive line-of-sight 
velocities on the higher- RA side, which is evident in Fig. |4] 

However, one should keep in mind that an observer seeing 
RX Boo from directions deviating considerably from the geo- 
centric direction would see a different set of maser components. 
Thus, Fig. Qis merely showing the three-dimensional positions 
of the maser components that we see from earth. Most of them 
probably would be too weak to be detected from the directions 
represented by Figs.|2](b) and (c). Of course a similar situation 
prevails for our present point of observation: there may be maser 
spots that we do not see because their emission may happen to 
be beamed in the wrong direction from our vantage point. 

Having adopted an expansion law it is of interest to investi- 
gate the associated time scale for the gas to travel through the en- 
velope. Thus, we have integrated the function 1 / v, where v is de- 
scribed by Eq. (O, along the radius r, avoiding r = ro = 1 .5 AU 
where v = 0. We find that it takes ~16.5 years for gas to travel 
through the zone (14-36 AU) where the H2O masers reside. 

6. Single dish data of SV Peg 

The SRV SV Peg was selected for monitoring because of the 
detection of several maser features in 1990 not seen in previous 
spectra. We suspected that this star may show different variabil- 
ity properties compared to RX Boo. The H2O maser was de- 
tected for the first time in 1974 by Dickinson ( |1976l l. showing a 
strong line of ~60 Jy at a line-of-sight velocity of 2 kms"'. The 
maser was re-detected by Engels et al. (1988) in 1984 show- 
ing two lines at 2 and 9 kms"' with similar flux densities of 
~12 Jy. The development of the maser emission since 1990 is 
shown in Fig. [8] The maser was observed also by Benson & 
Litde-Marenin (|T996] l in 1990 and by Takaba et al. ( fT994l ) in 
1991 yielding line profiles that were similar to ours. 

Our monitoring of SV Peg covered two periods 1990 - 1995 
and 2000 - 2007. Sample spectra are shown in Fig.|9] Until 1992 
the integrated maser emission was dominated by different single 
maser features with flux densities >10 Jy. Such features were 
not observed after 1992 and the integrated maser emission cor- 
respondingly declined abruptly during the first period. Also af- 
ter the monitoring restart in the year 2000 the maser emission 
remained weak or absent until mid 2006. All spectra taken are 
shown in the Appendix (Figs. lAiTj - IA.lOl l. 

The H2O masers of SV Peg appeared non-simultaneously 
over a velocity interval of to 16 km s"' . Actually, maser emis- 
sion in the outer parts of this interval (vios < 2 and >10 kms"') 
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was detected only during 1990 - 1993. The velocity range of the 
maser in general would be severely underestimated, if its deter- 
mination were based on individual spectra alone. 

Compared to RX Boo the spectral profiles of SV Peg are 
different. While for RX Boo generally maser emission was de- 
tected over the full velocity range and the individual features are 
strongly blended, the SV Peg maser features are usually well 
separated. This is in part due to the weakness of the SV Peg 
features, of which mostly only the brightest were detected. But 
even during the bright phases around 1991 and near the end of 
2006 most features are well separated suggesting that the num- 
ber of individual maser features is lower in S V Peg compared to 
RXBoo. 

Like we did for RX Boo (cf. Sect. l3.3b we made a line profile 
analysis for SV Peg by fitting multiple Gaussian components to 
the maser profiles. As before, we assumed that maser features in 
adjacent spectra with velocity differences ^0.5 kms"' belong to 
a unique physical maser component in the shell of SV Peg. All 
features were assigned to 9 maser spectral components labelled 
A, B,..., H. These components are listed in Tables IA3] and IA!41 
in the Appendix. For a description of the tables we refer to Sect. 
[331 

Except for features B and C we find that all spectral compo- 
nents are fairly well separated in velocity. Therefore the central 
velocities show in general a small scatter (^0.5 kms"'). Larger 
drifts in velocity (^1.5 kms"') are seen in component B indi- 
cating the presence of blends. Almost all features dominated 
the spectral profile at some epochs. The typical timescales for 
major intensity variations of individual components were 6 - 
12 months. A particularly well defined example for the tempo- 
rary appearance of a feature is component F (v = 11.3 kms"' 
in Fig. H Fig. |9l Date: 1991 March 30), which was detected 
for 1 1 months only. Noticeable is also the spectral feature H at 
-1-16.4 kms which was observed only twice in the beginning 
of 1992 (cf. FiglH Date: 1992 February 29). We consider this 
feature as real, although it was detected at the 2 - 3cr level only. 

Summarizing, the case of SV Peg shows that the H2O maser 
is composed of at least 8 components, where each of them dom- 
inated the profile for several months. At any given epoch, maser 
emission was only seen for a limited fraction (a few kms"') 
of the whole velocity span liable to exhibit a maser. Thus the 
non-detection of H2O masers in many SRVs might be merely 
a temporary effect, and the detection rate of 20% obtained by 
Szymczak & Engels (1995) among SRVs using two observing 
epochs separated by four months might be raised considerably 
with observations separated by years. 

7. Radio and optical variability in SRVs 

It has been reported that the variability of the H2O maser emis- 
sion of RX Boo is weakly correlated with the optical lightcurve, 
and that different maser features react to optical variations with 
different time delays (Hedden et al. |1991| ). The optical and the ra- 
dio lightcurve during the period 1987 - 2007 is shown in Fig.fTOl 
The optical data given as visual magnitudes were obtained from 
the American Association of Variable Star Observers (AAVSCQ; 
Hedden 2007). We rebinned the data into three-day intervals and 
interpolated over intervals without observations. The integrated 
radio fluxes were taken from Table lA.ll in the Appendix. 

A cross correlation between the optical and the radio 
lightcurves has been performed with delays of up to 2000 days 
(~6 years) following the suggestion of Lekht et al. ( 1200 11 1 that 

^ http://www.aavso.org/ 



the maser emission may respond to the passage of shock waves 
formed during a previous bright phase of the star. No signifi- 
cant signal has been found for any delay. However, the course of 
the lightcurves in Fig. [TO]gives the impression that a weak cor- 
relation without significant delay is present with a timescale of 
several years. In particular, both lightcurves show a minimum at 
2000 days (1992 - 1993). 

We searched also for particularly bright phases in the optical 
lightcurve before 1995, which could have caused the strong radio 
flare in 1995 March (f ^ 3000). Such opticafly bright phases 
occurred irregularly every few years, while the strength of the 
flare was unique during our observing period. It was therefore 
not possible to link the radio flare with a particular part of the 
optical lightcurve. 

We made a similar analysis of the correlation between opti- 
cal and maser variations for S V Peg with negative result. Lekht et 
al. (il999i ) observed the maser variations of another SRV RT Vir 
over the period 1985 - 1998 and found large amplitude varia- 
tions including two flares with flux densities > 1000 Jy. Although 
they do not discuss the optical behaviour of the star during this 
time, available AAVSO optical data do not show any magni- 
tude variations larger than the usual +0.5 mag around the mean 
Wvis ~ 8.9. We conclude therefore that there is basically no 
correlation between H2O maser radio data and optical data for 
SRVs. 

8. Discussion 

The long-term monitoring programme of the two SRVs RX Boo 
and SV Peg allows us to draw a few general conclusions on the 
variability properties of their H2O masers. 

- Maser emission occurs over a fixed velocity range Av, but 
not necessarily at all velocities simultaneously. There is no 
evidence that the velocity range varied systematically with 
time. Occasionally, maser emission outside the regular ve- 
locity interval may appear. 

- The individual maser features vary irregularly and any cor- 
relation with the optical variability of the stars is negligi- 
ble compared to these irregular variations. The intensity of a 
particular feature characterises this feature alone and is not 
representative for the maser emission as a whole. The inte- 
grated intensity obtained over the velocity range Ay is a more 
suitable measure to describe the maser "activity" at a given 
epoch. 

- The lifetime of individual maser components is of the order 
1 year, much less (5 - 10%) than the time needed for a mass 
element to cross the H2O maser shell. 

- The distribution of the maser spatial components in the cir- 
cumstellar shell of RX Boo has been found to be remarkably 
stable over the range of 6 years (1990 - 1995). The distribu- 
tion showed a persistent asymmetry, with emission coming 
always from the NE part and never from the SW part of the 
shell. The likelihood to find a new maser component in the 
shell, is higher in those parts, where a maser has been ob- 
served previously. 

8.1. RX Boo velocity range 

During the monitor programme the H2O maser velocity range of 
RX Boo never exceeded the interval -5 ^ vios ^ +5 km s ' (Fig. 
|2]i. This is well within the velocity range defined by CO emission 
extending from -8.3 to -1-10.3 km s ' according to Olofsson et al. 
ilOOlj) or -6.5 to +8.5 km s"' according to Teyssier et al. ( I2006I I, 
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Fig. 9. Sample H2O maser spectra of SV Peg showing typical profiles starting in 1990 and ending in early 2007. Note the varying 
flux density scale. 
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although the H2O interval is not symmetric with respect to the 
stellar velocity of v» = +1 kms"' as defined by the CO lines 
(see also Sect. I8.7l i. Larger velocity ranges have been observed 
in the past (Engels et al. |1993l l. Before 1983 the strongest feature 
was always detected at vios ~ 6 kms"' (see the compilation in 
Engels et al. 119881 1. which is well beyond the red edge of the cur- 
rent velocity interval. On one occasion, in 1982.2, emission was 
detected even over an interval - 11 < vios < -1-9 km s"' (Johnston 
et al. |1985l de Vegt et al. |19871 l. almost twice the range observed 
in 1980 (Dickinson & Dinger [T982j and after 1990. 

Except for this one occasion, the larger maser velocity range 
can be explained within our model (cf. Sect. |5]), as line-of-sight 
velocities |vios - v,| up to the adopted final outflow velocity 
Vf - 8.4 kms ' are permitted. According to Eq. ^ a larger 
I Vios - v,\ would require a larger outflow velocity v at the maser 
site and/or a larger 6. Following Eq. ^ the latter implies a 
smaller projected distance rp of the maser site. Because the out- 
flow velocity increases outwards, larger |vios - v,| are achieved 
at larger radial distances from the star and RX Boo may have at 
some time sustained H2O maser emission in the outer reaches of 
the shell. On the other hand, in the case of a smaller rp, maser 
components may have been present within the inner projected 
ring drawn in Fig. [H in which no emission was detected after 
1983. 

Maser emission at larger radial distances might have been 
possible due to a temporal increase in luminosity of the star 
allowing the rise of the temperature in the shell or due to the 
passage of a subshell with increased density. In the former case 
the minimum temperature required to excite H2O masers will be 
achieved at larger radii than normally. However, the mean opti- 
cal brightness of the star was not brighter before 1983 than af- 
terwards, making this explanation unlikely. In the latter case the 
size of the shell having the minimum density required to excite 
H2O masers might have been increased for some time. 

On the other hand, maser sites with smaller rp could have 
been triggered by local enhancements of the gas density or by 
temporary increases of gain lengths. Such sites could be located 
at any distance from the star, but we would be seeing a maser 
that has become more 'radial' . All these possible causes can be 
probed during future periods, if RX Boo should again display 
an extended velocity range, by searching for an increase of the 
outer shell radius, or by searching for maser components within 
the inner ring. 

At odds with our hypothesis is the 1982.2 observation, where 
the blueshifted emission exceeded the blue edge of the CO ve- 
locity range by 3 - 4 kms"'. This observation apparently has 
recorded a faster than usual outflow, crossing the H2O maser 
sheU by 1982. 

8.2. SV Peg velocity range 

The observed velocity interval in SV Peg is to +16.5 kms"'. 
This is to be compared to the CO emission velocity interval -2.5 
to -1-12.5 kms"' (Olofsson et al. 2002), which implies a stellar 
radial velocity v, = +5 kms"' and a final outflow velocity of 
7.5 km s"' . H2O maser emission at the most blue-shifted veloc- 
ities v'los < kms"' was never observed, while emission at the 
redshifted final outflow velocity was seen in 1990 - 1992 (spec- 
tral components F and G). This implies that the final outflow 
velocity is already reached within the H2O maser shell. 

As in the case of RX Boo before, also SV Peg showed 
emission outside the CO velocity range. Spectral feature H at 
H-16.4 kms"' (outflow velocity 11.4 kms"'), which was seen 
only in 1992, might have traced an outflowing mass segment 



with temporarily higher speed than is normally achieved in the 
CSE of SVPeg. 

8.3. H2O maser shell size of RX Boo 

Unless there are 100% tangentially beaming maser clouds at the 
outer periphery of the shell, the observed diameter 0p of the H2O 
maser shell will be smaller than the true diameter ^ because of 
projection eff'ects. As the most distant maser regions are usually 
weak, the observed diameter is certainly depending also on the 
sensitivity of the observations. Finally the diameter will depend 
on the assumptions made with respect to the location of the star 
relative to the observed distribution of maser clouds. 

For 1982 Johnston et al. (1985 ) give a projected diameter 
<pp = 200 mas compared to cpp = 140 mas in 1983 (Lane et al. 
iTWTT l and (pp = 150 mas in 1985 (Bowers et al. lT993l l. This ap- 
parently confirms the case, that the large velocity range in 1982 
was due (at least in part) to a more extended H2O maser shell. 
However, Colomer et al. ( i2000D determined (pp - 230 mas for 
1990 June and we fitted (pp - 180 mas for the bulk of the emis- 
sion features observed in 1990 - 1992, however, omitting outly- 
ing components at radial distances 100 < rp < 160 mas for the 
adopted stellar position. Usually the stellar position is adopted 
in the centre of the maser spot distribution, which leads to the 
smallest projected diameters possible. Our proposal for the po- 
sition of the star relative to the masers is offset by ~50 mas from 
the centroid of the maser emission, which leads to an unprece- 
dented high (pp - 320 mas. Thus the implied projected shell 
diameter is 5 1 AU. According to our model the component with 
the highest radial distance from the star has r ~ 36 AU, which 
implies that the true angular diameter of the shell might be as 
high as ip - 440 mas. So far, a variation of the H2O maser shell 
size of RX Boo cannot be determined, because the differences in 
the diameters derived by different authors are still dominated by 
data quality and model assumptions. 

8.4. The connection of maser variations with stellar 
variations 

Long-term monitoring of H2O maser emission has been made 
by Lekht et al. ( 1200 II and references therein) for various SRVs, 
Mira variables, and RSGs. For the objects with large-amplitude 
variations they find that the integrated H2O maser emission is 
correlated with the optical emission of the star The radio emis- 
sion lags behind the optical one by a few tenths of the period. 
They generally favour a model, in which the variations of the 
H2O masers are due to passages of shock waves generated on 
the stellar surface. The long travelling times (>10 years), how- 
ever, make it difficult to find correlations between the optical and 
radio lightcurves over such delays. 

However, the luminosity variations of the star certainly affect 
the excitation conditions in the whole maser shell. The response 
of H2O masers to the light variations in general is well docu- 
mented for Mira variables (e.g. Lekht et al. 120011 Brand et al. 
2002) and OH/IR stars (Engels et al. 19861 For the SRVs stud- 
ied here the variations of the integrated radio flux are largely 
caused by brightness variations of single maser features, and not 
by variations at all velocities simultaneously. We argue therefore 
that the H2O maser variations in SRVs are dominated by local ef- 
fects, such as variations of the velocity coherence length due to 
turbulence. The influence of the stellar luminosity variations on 
the maser variation is negligible in these stars. 
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8.5. The lifetime of individual maser clouds 

The changes seen in the VLA maps as well as the velocity shifts 
seen over time for spectral components determine a typical life- 
time of maser clouds in the range 0.5 - 3 years (Sect. |43] |. This 
is in excellent agreement with the estimates derived from H2O 
maser cloud imaging of four AGB stars of Bains et al. (120031 ). 
who argue that instability will destroy the clouds on such time 
scales. They derived sound crossing times of 1.5 - 3 years for 
the observed sizes of the clouds. 

Based on observations with the VLBA Imai et al. ( 120031) de- 
rived for RT Vir clouds considerably shorter lifetimes of 1 - 2 
months. This order of magnitude shorter lifetime compared to 
ours is probably due to different spatial resolutions. The number 
of clouds identified in a particular maser shell is significantly 
increased when the spatial resolution is increased. For exam- 
ple. Marvel ( 11996 !) distinguished >50 maser clouds in VLBA 
RX Boo H2O maser maps, while we were restricted to <15. 
Thus, the clouds identified here for RX Boo are certainly com- 
posed of several sub-components, where each of them might sur- 
vive on the timescale of months only. Our larger lifetimes ap- 
ply then to spatially closely located groups of maser clouds (our 
maser spatial components), of which the spectral components in 
single-dish spectra are made. The dependence of the lifetimes 
on the size of the maser clouds is also evident from observations 
of RSGs, where the clouds have physical diameters a magnitude 
larger than are observed in AGB stars, and lifetimes of the order 
of 10 years (Richards et al. 1998 and 1999). 

The typical gas crossing time through the H2O maser sup- 
porting shell in RX Boo (thickness ~22 AU) is ~16.5 years 
(Sect.|5]l. The maser clouds most probably do not survive a full 
crossing of this shell. Due to the limited lifetimes of the maser 
clouds proper motion measurements require therefore an observ- 
ing frequency high enough to sample a given cloud at least twice. 
At the highest resolutions achievable today, recurrent observa- 
tions with intervals less than a month are reasonable (see Imai et 
al. [2003]) . 

8.6. The beaming directions of masers in the RX Boo shell 

The distribution of maser components in Fig. H] tells us that the 
individual masers are radiating non-isotropically in the sense 
that their maximum output is in a direction that does not deviate 
very much from the tangential direction. Our model of the shell 
expansion has given us values of the angle 6 (the aspect angle) 
for every maser component. The distribution function of these 
values has a maximum at zero degrees and an extent from about 
-60° to -1-60°. In other words, the masers are not all in a plane 
perpendicular to the line of sight [see Figs.|2lb) and (c))]. 

Now, imagine for a moment that all masers were radiating 
in the tangential direction and that we therefore see many of 
them from an angle that deviates from the direction of max- 
imum amplification. Furthermore, assume that all the masers 
have the same intrinsic strength. Then by plotting the observed 
flux density against the aspect angle we would get the 'beam 
shape' of the 'unit maser'. Of course there is probably a spread 
in the direction of maximum amplification from the tangential 
direction. Also there is certainly a considerable spread in the in- 
trinsic strength of the masers. These spreads would destroy to a 
large extent the relation between flux density and aspect angle. 
However, it is conceivable that a weak correlation would remain. 

In Fig. [TT] the observed flux density is plotted against the 
absolute value of the aspect angle. It is obvious that the highest 
fluxes are found among the smallest aspect angles. If the two 



data points with the very high fluxes close to 20° are disregarded 
the highest flux indeed is found close to 0°. In addition, the upper 
envelope of the scatter plot is decreasing from a maximum close 
to 0° to a minimum close to 60°. The fact that this correlation is 
disturbed by the very strong masers at 20° must be explained by 
the disparate intrinsic maser strengths. 

Fig.[TT]is an additional evidence, albeit a weak one, for the 
H2O masers in the CSE of RX Boo to be predominantly tangen- 
tially beamed. However, the typical beaming angle is impossible 
to extract from the data as we cannot disentangle the deviation of 
the line of sight from the direction of the maximum amplification 
and the deviation of the direction of the maximum amplification 
from the tangential direction. 

8. 7. Maser shell as)^mmetry 

So far, the H2O maser shell has been assumed to be radially sym- 
metric. However, there is a prominent asymmetry present in the 
distribution of H2O maser emission of RX Boo, which mani- 
fests itself during 1990 - 1995 by the absence of emission in the 
SW part of the shell and by the NW-SE velocity gradient (Sect. 
14.41 ). How persistent is this asymmetry? Interferometric observa- 
tions of the H2O masers of RX Boo have been made a number 
of times, but the first sensitive map available is from VLA ob- 
servations in 1985 (Bowers et al. 1993). An E-W velocity gra- 
dient is evident also in this map, but in addition there is some 
emission in the SW part of the shell. VLBA observations taken 
in 1995 January and July are available from Marvel ( 1996 ), but 
they are difficult to compare with our VLA maps because all the 
blueshifted emission with vios < +0.5 kms"' was not detected. 
This emission seems to have been resolved out. Nevertheless, we 
conclude that the incomplete ring-like structure and the velocity 
gradient was present at least over eleven years between 1985 and 
1995, which is about 70% of the gas crossing time in the H2O 
maser-supporting shell of RX Boo. 

Evidence for asymmetry is coming also from the H2O- 
maser velocity interval, which is centred on Vc = kms"' 
(Fig. |2|i, while the CO-emission velocity interval is centred on 
v, = -Hi kms"' (Teyssier et al. 2006). The blueshifted projected 
expansion velocities extend to higher values |vios-v»| * 6kms"' 
than the redshifted velocities |vios - v,| 4 km s"' . This velocity 
asymmetry has persisted over the full monitoring time range of 
20 years. 

In Fig. [12] the upper envelopes of the spectral profiles of 
RX Boo and SV Peg are shown. These envelopes show the 
highest flux densities measured during the monitoring period at 
each velocity. Contrary to SV Peg, where the highest flux den- 
sities are more or less on a similar level across the maser ve- 
locity range, the profile of RX Boo is dominated by emission at 
~2.5 kms"' (spectral component I). Emission close to this ve- 
locity was present over the full monitoring period (see Table |A.4l 
in the Appendix), emphasizing the longevity of the asymmetries. 

Asymmetries present in individual observations can easily 
be explained as due to shells unevenly filled with maser clouds. 
The clouds themselves are explained by Richards et al. ( 119991 ) 
as clumps of enhanced density, which form as condensations in 
the extended stellar atmosphere. Due to the randomness of the 
formation process of such clumps, parts of the shell might not 
contain clumps at a given epoch. But we have observed the per- 
sistence of the asymmetric emission structure and the velocity 
gradient over a time range more than half of the crossing time 
through the shell. And, the asymmetry of the H2O maser veloc- 
ity interval outlasted the 16.5 years crossing time, showing that 
the morphology of the maser site distribution is not due to ran- 
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dom fluctuations. It is therefore more likely that the persistent 
asymmetries are evidence for a non spherically-symmetric mass 
loss process in RX Boo. The asymmetric distribution of maser 
emission components found in RX Boo is not untypical. ^From 
VLA-observations of H2O masers in 13 Mira variables and 3 
SRVs, Bowers & Johnston ( 119941 ) concluded that the emission 
is best modelled by "clumps, segments, or filaments in an asym- 
metric but ellipsoidal-type geometry". 

For a strong maser to occur, some very special (and rare) 
physical conditions have to be fulfilled: the pumping of the mas- 
ing transition has to win over thermal de-excitation processes 
and the coherent amplification path length has to be sufficiently 
long for the maser flux density to be detected. Therefore, an 
alternative explanation of the asymmetric distribution of H2O 
masers in the envelope of RX Boo might be provided by a non- 
spherical distribution of maser conditions. 

A possible explanation of such a distribution might be 
found in the photospheric structure of the star Giant and su- 
pergiant stars have extremely deep convective layers and as 
a result the convection cells of the photosphere are expected 
to be gigantic (cf. Schwarzschild 1975). Freytag ( I2003I I has 
made a 3 -dimensional radiative hydrodynamic model of the 
SRV Betelgeuse {a Orionis; spectral class M21ab) that simu- 
lates high-resolution observations of surface features (cf. Young 
et al. '2000) and typical irregular flux variations of this star (cf. 
Goldberg 1984 ). However, we are unaware of any published sim- 
ulations of pulsating giants or of AGB stars of more moderate 
sizes. Nevertheless, it is conceivable that the random character 
of the convection elements is responsible for irregular light vari- 
ations or for an irregular modulation of periodic variations. It 
is conceivable too that the convection elements to a certain de- 
gree are preserved as some kind of 'blobs' in the stellar wind. 
As these blobs reach the H2O maser region, a larger blob prob- 
ably creates longer amplification path lengths than a number of 
smaller blobs. In addition such blobs will expand considerably 
during their travel outwards, not only due to the divergence but 
also radially due to the acceleration. In the case of RX Boo a 
convective blob is expected to increase its volume by a factor of 
~50 during its journey from the photosphere to the inner radius 
of the H2O maser region. 

However, the real mystery is the longevity of the asymme- 
try of the H2O maser region of RX Boo. In the convection-wind 
picture this implies that the lifetime of convection cells either be 
long or that the occurence of large cells persists over long time 
periods. In the Betelgeuse simulation by Freytag (I2003I I the life- 
time of convection cells seems to be of the order of one year 
and the likelihood of occurence of a cell seems to be constant 
over both time and stellar surface. If the typical lifetime of con- 
vection cells in the RX Boo photosphere is of the same order, it 
would take about one year for the corresponding blob to getting 
detached from the photosphere and becoming part of the stellar 
wind. However, when this blob reaches the H2O maser region 
(~14 AU from the star), the inflation due to acceleration and di- 
vergence would increase the time for the blob to pass a fixed 
point by a factor of about 5. 

9. Conclusions 

Water vapour maser emission from CSEs is often strong enough 
to be observed with small radio telescopes even in rather dis- 
tant late-type stars. Because of its variability, however, the in- 
formation to be derived from a single maser spectrum is limited. 
Neither the intensities nor the velocities are necessarily repre- 
sentative. To understand the cause of the strong variability of the 



masers in SRVs we have monitored RX Boo and SV Peg, two 
stars representative of their class, with single-dish radio tele- 
scopes over a period of two decades. The interpretation of the 
variations seen in single-dish spectra is usually compromised by 
blending of the maser peaks due to the overlap of two or more 
maser clouds in velocity space. To overcome this limitation, we 
obtained in addition for RX Boo four interferometric observa- 
tions taken over a time interval of almost three years. 

We find that the H2O masers in RX Boo are located in an 
incomplete ring of thickness 22 AU, which is filled only in part. 
The time that the gas needs to cross this ring is 16.5 years. This 
ring-like distribution of maser spots remained stable over at least 
1 1 years. Compared to these timescales, the survival time of indi- 
vidual maser clouds is much shorter having a typical duration of 
~1 year The variability of H2O masers is then due to the emer- 
gence and disappearance of shortlived maser clouds. The vari- 
ations in intensities of these clouds are not dependent on varia- 
tions of the stellar luminosity, as we find no correlation between 
the optical and the radio lightcurve. 

It follows from the interpretation of single-dish spectra that 
the strength of individual maser peaks or the details in the shape 
of the profile are snap-shots reflecting the distribution of the 
maser clouds and their luminosities at the time of observation. 
If RX Boo and SV Peg are typical SRVs, the lifetime of these 
clouds in SRVs is limited to about 1 year. Also stellar radial ve- 
locities and shell outflow velocities will be biased, if the maser 
cloud distribution is not uniform in the shell. 

For RX Boo, the persistent overall distribution of the maser 
emission regions indicates that the underlying mass outflow 
from the star or the distribution of maser conditions are not 
spherically symmetric. Whether deviations from spherical sym- 
metry are present in the H2O maser regions of SRVs in general, 
will have to be answered by new monitoring programmes rely- 
ing on interferometry for a larger sample of stars. 
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Table 3. Spatial and spectral components of RX Boo 
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Fig. 7. The three-dimensional distribution of maser components 
and their associated velocities according to the expansion model 
given by Eq. ((S). The distribution is seen from three cardinal 
directions: on the sky (a), 'from the side' (b) and 'from above' 
(c). The maser components are represented by the same symbols 
as earlier. In addition the projected expansion velocities and the 
line-of-sight velocities are shown in the same scale in units of 
km s ' as the distance scale in units of AU. 
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Fig. 8. Flux density of the H2O maser of S V Peg as a function of 
velocity and time. The time axis covers the years 1990 - 1995 
(bottom) and 2000 - 2007 (top). The time stamp f = corre- 
sponds to 1990 Februai-y 17 (J.D. 2447933). The dotted hori- 
zontal lines mark the actual observations averaged in 4-day in- 
tervals. Flux densities in between the observations were inter- 
polated linearly. The dashed vertical line marks the stellar sys- 
temic velocity (5 kms"') as given by Kerschbaum & Olofsson 
( I1999I I. The spectra were resampled to a resolution of 0.3 km s ' . 
Contours in the lower panel are drawn at levels 3.00, 4.19, 5.86, 
8.18, 11.43, 15.97 Jy (logarithmic scale) and in the upper panel 
at 3.00, 4.40, 6.45, 9.45, 13.86, 20.32 Jy. These levels are indi- 
cated by black and white lines in the flux scales along the right 
side of the two maps. (For a colour version of this figure, see the 
electronic edition.) 
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radio flux has been offset by an arbitrary amount. 
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Fig. 11. The observed flux density as a function of the absolute 
value of the aspect angle (6). The epochs are represented by the 
same symbols as in Figs. |4] and |5] 
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Fig. 12. Upper envelope of all spectral profiles of RX Boo (top) 
and SV Peg (bottom). In RX Boo the spectra containing the 
flare (1995 January - May) were omitted. The flare occured at 
~3.5 kms"' close to the strongest peak in the upper envelope 
and reached a peak flux of 1600 Jy. The upper envelope is de- 
fined as the highest flux density ever achieved in velocity bins of 
0.3 km s ' during the monitoring program. 
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Table A.l. Maser components of RX Boo. The columns give the observing date in Gregorian and Julian (JD) format, the noise of the spectra (rms) 
the integrated flux 5i in units of 10"^- W nr^ 
A - E in Tabled 



and the velocities Vp and flux densities 5p of maser peaks associated with the spectral components 
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Table A.l. continued. 
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-6.0 - 0.0 « 


2-4 






















20.06.02 


2452446 


fl 1 ^ 

U. i J 


^19 

J iZ 


— 


— 






4.0 


-0.0 
























26.06.02 


2452452 


1 1 1 


A'XA 


— 


— 


— 


— 






— 


— 


— 


01.10.02 


2452549 


fl Q^> 


48^ 
^o J 


— 


— 


— 


— 






— 


— 


— 


24.10.02 


2452572 


fl QQ 
U.yy 


987 
Zo / 


— 


— 


— 


— 






— 


— 


— 


13.12.02 


2452622 


1 An 
1.4/ 


^11 

J 1 1 


— 


— 


— 


— 






— 


— 


— 


19.12.02 


2452628 


1 .J J 


48fl 
4oU 


— 


— 


— 


— 






— 


— 


— 


15.01.03 


2452655 


i .J J 


488 
4oo 


— 


— 


— 


— 






— 


— 


— 


02.04.03 


2452732 


1 fl^ 
i .UJ 


49^ 
4ZJ 


— 


— 


— 


— 






— 


— 


— 


18.11.03 


2452962 


fl QA 
U.oO 


7^7 
/J / 


— 


— 


— 


— 






— 


— 


— 


28.01.04 


2453033 


fl Sfl 
U.oU 


41 
^lO 


-5.1 


3.5: 


— 


— 






— 


-1.4 


22.9 


31.03.04 


2453096 


fl 8^ 
W.o J 


jyj 


— 


— 


— 


— 






— 


-1.3 


21.1 


11.05.04 


2453137 


fl A'7 
U.u/ 


Z40 


— 


— 


— 


— 






— 


-1.3 


4.6 


18.06.04 


2453175 


fl ^/l 
U.j4 


9 1 A 
ZiO 


— 


— 


— 


— 






-2.0 - 0.0 


a 2 




17.09.04 


2453266 


1 flfl 
1 .uu 


^A9 
JOZ 


— 


— 


— 


— 






— 


— 


— 


18.12.04 


2453358 


114 
1.1^ 


jZj 


— 


— 


-3.6 


7.3 






-2.0 - 0.0 


X 2 




12.01.05 


2453383 


fl Q^^ 


447 


— 


— 


-3.6 


8.3 






-2.0 - 0.0 


K 2 




15.02.05 


2453417 


fl 7/1 
U. /4 


^/l 8 
j4o 


— 


— 


-3.9 


3.5 


-2.8 2.9 




— 


-0.4 


3.2 


13.04.05 


2453474 


1 01 


7Sfl 
/ JU 


— 


— 


— 


— 







— 


— 


— 


20.06.05 


2453542 


1.56 


394 


— 


— 


— 


— 






— 


— 


— 


11.07.05 


2453563 


1.13 


375 


— 


— 


— 


— 






— 


— 


— 


22.11.05 


2453697 


1.26 


703 


— 


— 


— 


— 


-2.6 4.0: 




— 


— 


— 


14.Uz.UO 




1.32 


940 




















07.04.06 


2453833 


1.22 


527 




















05.07.06 


2453922 


1.01 


380 




















01.09.06 


2453980 


1.33 


436 




















17.10.06 


2454026 


0.62 


436 




















01.12.06 


2454071 


1.01 


542 




















17.01.07 


2454118 


0.94 


320 
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Table A.2. Maser components of RX Boo. The columns give the observing date in Gregorian and Julian (JD) format, and the velocities Vp and flux 
densities 5p of maser peaks associated with the spectral components F ~ K in Table[3] 



G 



Date 



JD 



30.03.87 
22.06.87 
05.09.87 
12.09.87 
05.11.87 
11.02.88 
19.02.88 
02.05.88 
13.07.88 
13.10.88 
31.01.89 
21.03.89 
04.12.89 
06.02.90 
17.02.90 
31.03.90 
17.04.90 
11.05.90 
15.07.90 
21.10.90 
18.01.91 
29.01.91 
30.03.91 
01.05.91 
17.05.91 
25.10.91 
18.01.92 
29.02.92 
18.04.92 
05.07.92 
01.09.92 
15.10.92 
24.10.92 
22.12.92 
26.01.93 
21.04.93 
13.05.93 
03.11.93 
30.11.93 
08.03.94 
14.04.94 
08.09.94 
28.10.94 
18.01.95 
09.03.95 
12.03.95 
13.03.95 
27.03.95 
09.04.95 
03.05.95 
03.06.95 
04.06.95 
05.06.95 
23.06.95 
18.09.95 
18.01.96 



H I 

[kms"', Jy] . 



J 



K 



2446885 
2446969 
2447044 
2447051 
2447105 
2447203 
2447211 
2447284 
2447356 
2447448 
2447558 
2447607 
2447865 
2447929 
2447940 
2447982 
2447999 
2448023 
2448088 
2448186 
2448275 
2448286 
2448346 
2448378 
2448394 
2448555 
2448640 
2448682 
2448731 
2448809 
2448867 
2448911 
2448920 
2448979 
2449014 
2449099 
2449121 
2449295 
2449322 
2449420 

2449457 

2449604 

2449654 

2449736 

2449786 ■ 

2449789 

2449790 ■ 

2449804 

2449817 

2449841 

2449872 

2449873 

2449874 

2449892 

2449982 

2450101 



0.1 8.5 



0.1 9.7 



0.1 6.2 



0.9 
1.4 



1.3 
0.6 

0.4 
0.5 



22.7 
22.5 
1.6 
- 1.6 



11.0 

11.9 
9.9 



1.4 


71.0 


2.4 


212.8 


1.3 


24: 


2.5 


173.6 


1.3 


27: 


2.5 


261.9 


1.5 


32: 


2.5 


294.0 


1.3 


36.3 


2.5 


172.4 


1.6 


41: 


2.6 


215.8 


1.6 


44: 


2.7 


232.4 


1.2 


39: 


2.8 


238.7 


1.7 


34: 


2.9 


108.0 


1.9 


61.3 


3.0 


136.8 


1.9 


35.6 


2.9 


93.2 


32.2 




2.9 


92.6 


9.7 




2.9 


34.5 


1.4 


35.3 


2.7 


89 


1.6 


35.8 


2.9 


64.7 


1.6 


30.7 


2.8 


52.1 


1.5 


39.4 


2.7 


48.9 


1.6 


46.3 


2.8 


47.3 


1.5 


13.2 


2.5 


28.7 


1.8 


35.4 


2.8 


16.4 




4.1 40.2 
4.0 18.6 
4.0 18: 



4.9 
5.1 
5.1 



5.6 
9.5 
4.3 



4.5 2.3 



2.8 
8.1 
5.1 
1.8 
1.9 



0.0 
0.0 



0.4 23.3 
0.6 10.5 
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H 



I 



K 



Date 



JD 



01.03.96 

03.09.96 

14.03.97 

06.05.97 

12.10.97 

21.10.97 

15.12.97 

29.01.98 

19.03.98 

10.04.98 

11.05.98 

02.02.99 

24.06.99 

27.10.00 

20.12.00 

07.05.01 

18.09.01 

24.10.01 

19.03.02 

18.04.02 

20.06.02 

26.06.02 

01.10.02 

24.10.02 

13.12.02 

19.12.02 

15.01.03 

02.04.03 

18.11.03 

28.01.04 

31.03.04 

11.05.04 

18.06.04 

17.09.04 

18.12.04 

12.01.05 

15.02.05 

13.04.05 

20.06.05 

11.07.05 

22.11.05 

14.02.06 

07.04.06 

05.07.06 

01.09.06 

17.10.06 

01.12.06 

17.01.07 



2450144 
2450330 
2450522 
2450575 
2450734 
2450743 
2450798 
2450843 
2450892 
2450914 
2450945 
2451212 
2451354 
2451845 
2451899 
2452037 
2452171 
2452207 
2452353 
2452383 
2452446 
2452452 
2452549 
2452572 
2452622 
2452628 
2452655 
2452732 
2452962 
2453033 
2453096 
2453137 
2453175 
2453266 
2453358 
2453383 
2453417 
2453474 
2453542 
2453563 
2453697 
2453781 
2453833 
2453922 
2453980 
2454026 
2454071 
2454118 



"p "p Jp 
_[kms-\Jy] 



0.0 
0.0 



0.9 15.2 
15.7 
4.4 



0.5 

0.2 



0.1 8.6 



2.3 



2.5 
2.6 
2.0 
1.9 
1.9 
2.0 
1.8 
1.6 



2.6 
2.5 



6.5 



9 1 
z. 1 


O.J 


1 


<o 


L.D 




1 1 
z.z 


L l.J 


Z.J 


J.O 


z.z 




1 c 

l.J 


J.J 


2.6 


8.8 


2.6 


7.9 


2.6 


10.7 


2.2 


17.4 


2.2 


27.1 


2.3 


26.7 


2.2 


31.5 


2.3 


12.2 


2.4 


11.8 


2.7 


24.3 


2.8 


14.3 


3.0 


22.8 


3.0 


15.8 



7.4 

5.9 
<4.5 

4.2 
15.8 

7.7 
11.7 
14.3 



2.6 

2.6 

34.0 

18.1 

2.7 



2.8 

2.9 

3.2 

3.2 

3.4 

3.4 

3.5 

3.5 

3.5 

3.3 

3.5 

3.2 

3.0 



3.0 
3.3 
3.6 
3.5 
3.5 
3.5 
3.5 
3.5 
3.7 
3.6 
3.4 
3.3 
3.4 
3.4 
3.5 
3.5 
3.6 
3.7 
3.6 
3.4 
3.3 
3.1 
3.2 
3.3 
38.3 
86.3 
3.1: 
3.4: 
20.1 
2.9 
2.9 
2.6 



13.8 
16.4 
23.0 
36.5 
17.7 
13.3 
17.2 
11.8 
12.0 
9.2 
14.9 
20.5 
51.4 



10.7 
10.3 
84.6 
89.3 
30.8 
25.5 
19.4 
13.3 
25.2 
25.8 
39.2 
40.9 
111.8 
17.6 
13.8 
12.0 

9.6 
11.8 

8.8 
13.4 
19.1 
31.3 
19.5 
16.5 



14.3: 
10.6: 



20.5 
27.5 
13.7 



4.1 



4.9 16.5 



<7 



4.4 41.7 

4.4 8.5 

4.0 8.4 

3.9 18.6 



4.5 15.5 

4.6 17.8 
4.6 18.4 
4.8 13.3 
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Table A.3. Maser components of S V Peg. The columns give the observing date in Gregorian and Julian (JD) format, the noise of the spectra (rms), 
the integrated flux 5i in units of 10"^- W irT^, and the velocities Vp and flux densities 5p of maser peaks associated with the spectral components 
A - D as described in Section[6] 



A B C D 





Date 


JD 


rms 

iJyi 


Si 


K 


S n ^'n 


Sr, 

[km s 


K 

, Jyj 


Sn V 


n S r) 




17.02.90 


2447940 


0.21 


142 


— 


— 


2.2 


22.9 


— 


— 


— 


— 


01.04.90 


2447983 


0.19 


82 


— 


- 


2.0 


3.2 


- 


— 


— 


- 


13.05.90 


2448025 


0.15 


95 


— 




2.2 


3.8 


- 


- 


- 


- 


25.10.90 


2448190 


0.20 


92 


0.5 


0.7 


2.3 


4.3 


— 


— 


4.5 


0.9 


19.01.91 


2448276 


1.42 


<50 


— 


— 


2.4 


3.4 


— 


— 


— 


— 


30.03.91 


2448346 


0.16 


122 


0.3 


0.4 


2.3 


1.7 


- 


— 


4.5 


0.9 


01.05.91 


2448378 


0.22 


136 


0.3 


0.5 


2.3 


1.8 


— 


— 


4.3 


1.4 


22.05.91 


2448399 


1.67 


145 


— 


- 


— 


— 


— 


— 


— 


— 


25.10.91 


2448555 


1.99 


219 


— 


- 


2.1 


5.3 


- 


— 


4.5 


7.3 


18.01.92 


2448640 


0.19 


89 


0.3 


0.8 


2.0 


1.9 


— 


— 


4.7 


3.9 


29.02.92 


2448682 


0.15 


89 


1.0 


0.3: 


2.1 


0.9 


- 


— 


4.7 


4.5 


19.04.92 


2448732 


1.01 


327 


- 


- 


— 


— 


- 


— 


4.4 


18.6 


02.09.92 


2448868 


0.20 


102 


0.5 


1.2 


— 


— 


3.6 


12.2 


— 


— 


16.10.92 


2448912 


0.89 


252 


0.6 


6.0 


— 


— 


3.4 


24.4 


— 


— 


22.12.92 


2448979 


0.20 


180 


0.3 


16.1 


1.7 


1.5 


3.5 


3.7 


4.9 


0.7: 


25.01.93 


2449013 


0.48 


101 


0.4 


11.7 


1.9 


2.0: 




4.2 


1.4 




20.04.93 


2449098 


0.13 


84 


0.3 


4.1 


1.5 


2.6 


3.5 


2.2 


4.8 


2.2 


02.11.93 


2449294 


0.65 


74 


0.7 


1.7 


- 


- 


3.3 


2.4 


4.8 


2.5 


29.11.93 


2449321 


0.33 


59 




1.2 


1.0 




3.0 


2.6 


4.7 


4.3 


17.01.94 


2449370 


0.47 


64 


— 


— 


_ 


_ 


_ 


_ 


4.5 


6.0 


08.03.94 


2449420 


0.22 


30 


- 


- 


1.6 


0.7 


3.3 


0.9 


4.5 


2.3 


15.04.94 


2449458 


1.16 


<50 


— 


— 


— 


— 


— 


— 


— 


— 


18.01.95 


2449736 


0.52 


<50 


— 


— 


2.0 


1.7: 


— 


— 


— 


— 


09.03.95 


2449786 


0.09 


19 


— 


— 


2.0 


1.0 




4.1 


0.4 




24.06.95 


2449893 


0.23 


17 


— 


— 


1.9 


0.7 




4.0 


1.0 




27.10.00 


2451845 


1.15 


<50 


— 


— 


_ 


_ 


_ 


_ 


_ 


_ 


20.12.00 


2451899 


0.94 


<50 


— 


— 


— 


— 


— 


— 


— 


— 


07.05.01 


2452037 


1.09 


<50 


— 


— 


— 


— 


— 


— 


— 


— 


18.09.01 


2452171 


0.87 


<50 


— 


- 


— 


— 


— 


— 


— 


— 


24.10.01 


2452207 


1.21 


<50 


— 


- 


- 


- 


- 


- 


- 


- 


20.03.02 


2452354 


1.42 


<50 


— 


- 


— 


— 


— 


— 


— 


— 


24.04.02 


2452385 


1.24 


65 


— 


— 


— 


— 


— 


— 


— 


— 


27.06.02 


2452453 


2.87 


<50 


— 


— 


— 


— 


— 


— 


— 


— 


01.10.02 


2452549 


0.94 


75 


- 


- 


2.8 


1.8: 


- 


— 


— 


- 


29.10.02 


2452577 


1.62 


<50 


— 


- 


— 


— 


— 


— 


— 


— 


19.12.02 


2452628 


1.11 


<50 


— 


- 


3.2 


1.7: 


- 


— 


- 


- 


14.01.03 


2452654 


1.05 


<50 


- 


- 


— 


- 


- 


— 


— 


— 


03.04.03 


2452733 


1.47 


<50 


- 


- 


— 


— 


— 


— 


— 


— 


19.11.03 


2452963 


1.37 


<50 


— 


— 


3.1 


2.1: 


— 


— 


— 


— 


28.01.04 


2453033 


0.76 


<50 


— 


— 


2.9 


1.5: 


— 


— 


— 


— 


02.04.04 


2453099 


0.80 


<50 


- 


- 


2.3 


1.9: 


— 


— 


— 


— 


14.05.04 


2453140 


0.99 


<50 


— 


— 


— 


— 


— 


— 


— 


— 


20.06.04 


2453177 


2.29 


<50 


— 


- 


— 


— 


— 


— 


— 


— 


17.09.04 


2453266 


1.88 


<50 


— 


— 


— 


— 


— 


— 


— 


— 


20.12.04 


2453360 


0.77 


54 


— 


- 


2.8 


2.0 


- 


— 


— 


- 


11.01.05 


2453382 


0.90 


<50 


— 


— 


_ 


_ 


_ 


_ 


_ 


_ 


12.03.05 


2453442 


1.08 


56 


















15.04.05 


2453476 


1.02 


104 






2.9 


3.0 










21.06.05 


2453543 


1.94 


<50 


















12.07.05 


2453564 


1.46 


81 


















22.11.05 


2453697 


1.33 


175 






2.6 


6.2 










15.02.06 


2453782 


1.61 


193 






2.4 


4.7 










09.04.06 


2453835 


2.57 


179 






2.5 


6.0 










07.07.06 


2453924 


2.49 


216 






2.1 


5.3 






6.1 


4.6 


01.09.06 


2453980 


1.00 


402 






2.2 


6.6 






5.7 


8.2 


17.10.06 


2454026 


0.71 


518 






2.1 


15.5 






5.7 


13.7 
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A B C D 





Date 


JD 


rms 


Vp 


Sp Vp 




Vp ^p Vp 








[Jy] 






[kms" 






30.11.06 


2454070 


1.47 


260 




2.0 


15.0 




5.5: 3.4: 


16.01.07 


2454117 


1.29 


322 




2.1 


35.8 
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Table A.4. Maser components of SV Peg. The columns give the observing date in Gregorian and Julian (JD) format, and the velocities Vp and flux 
densities 5p of maser peaks associated with the spectral components E - H as described in Section[6] 





Date 


JD 


E 

Vp 5p Vp 


F 

5p 
[kms" 


G 
Jyl 


5p Vp 


H 

S p 












17.02.90 


2447940 


7.7 






- 


12.9 


0.9 


- 


01.04.90 


2447983 


7.7 


1 8 
i . o 




- 


12.9 


2.3 


- 


13.05.90 


2448025 


- 






- 


- 


- 


- 


25.10.90 


2448190 


7.7 






- 


13.4 


0.4: 


- 


19.01.91 


2448276 


7.7 


d A 




- 


- 


- 


- 


30.03.91 


2448346 


7.4 


1 f, 

i .u 


1 1 A 


12.7 


13.4 


0.9 


- 


01.05.91 


2448378 


7.3 


1 ? 


1 1 A 


14.2 


13.4 


1.0 


- 


22.05.91 


2448399 


- 






11.8 


- 


- 


- 


25.10.91 


2448555 


- 




1 1 7 


5.7 


13.7 


2.6 


- 


18.01.92 


2448640 


- 






0.5 


13.7 


0.5 


16.4 


29.02.92 


2448682 


- 






- 


13.8 


0.5 


16.6 


19.04.92 


2448732 


- 






- 


- 


- 


- 


02.09.92 


2448868 


7.5 


1 7 




- 


- 


- 


- 


16.10.92 


2448912 


7.5 


7 7 




- 


- 


- 


- 


22.12.92 


2448979 


7.3 






- 


- 


- 


- 


25.01.93 


2449013 


- 
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Table A.4. continued. 
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E F G H 
Date JD Vp 5p Vp 5p Vp Sp Vp Sp 
[kms"', Jy] 

30.11.06 2454070 8.6 4.1 

16.01.07 2454117 _______ 



List of Objects 

'RX Boo' on page 13 
'SVPeg' on pagei 
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Fig. A.l. H2O maser spectra of RX Boo. Note the varying flux density scale. 
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Fig. A.2. H2O maser spectra of RX Boo (continued). Note the varying flux density scale. 
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Fig. A.3. H2O maser spectra of RX Boo (continued). Note the varying flux density scale. 
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. A.4. H2O maser spectra of RX Boo (continued). Note the varying flux density scale. 
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Fig. A.5. H2O maser spectra of RX Boo (continued). Note the varying flux density scale. 
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Fig. A.6. H2O maser spectra of RX Boo (continued). Note the varying flux density scale. 
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Fig. A.7. H2O maser spectra of SV Peg. Note the varying flux density scale. 
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Fig. A.8. H2O maser spectra of SV Peg (continued). Note the varying flux density scale. 
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Fig. A.9. H2O maser spectra of SV Peg (continued). Note the varying flux density scale. 
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Fig. A.IO. H2O maser spectra of SV Peg (continued). Note the varying flux density scale. 



